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Abstract: The presence of two light higgsinos nearly degenerate in mass is one of the 
important characteristics of supersymmetric models meeting the naturalness criteria. Prob¬ 
ing such higgsinos at the LHC is very challenging, in particular when the mass-splitting 
between them is less than 5GeV. In this study, we analyze such a degenerate higgsino 
scenario by exploiting the high collinearity between the two muons which originate from 
the decay of the heavier higgsino into the lighter one and which are accompanied by a 
high-pT QCD jet. Using our method, we can achieve a statistical signihcance ^ 2.9 cr as 
well as S/B ^ 17% with an integrated luminosity of 3000 fb“^ at the 14TeV LHC, for the 
pair production of higgsinos with masses 124 GeV and 120 GeV. A good sensitivity can be 
achieved even for a smaller mass-splitting when the higgsinos are lighter. 
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1 Introduction 

One of the key theoretical motivations for low-energy supersymmetry (SUSY) is that it 
provides a framework in which a light Higgs boson can be obtained without invoking unnat¬ 
ural fine-tuning of theory parameters. However, the Higgs boson discovered recently [1, 2] 
at the Large Hadron Collider (LHC) has a mass around 125 GeV and signal rates consistent 
with those predicted by the Standard Model (SM). These properties of the Higgs boson, 
in conjunction with the non-observation of supersymmetric particles, have resulted in ex¬ 
cluding large portions of the parameter space of the Minimal Supersymmetric Standard 
Model (MSSM) where the naturalness criteria are satisfied. If the observed Higgs reso¬ 
nance is to be identified with the lightest CP-even Higgs boson, /i, of the MSSM, TeV-scale 
SUSY-breaking masses and/or multi-TeV soft trilinear coupling parameters are necessary, 
so that the Higgs boson mass can be enhanced sufficiently via radiative corrections [3-8]. 
Furthermore, null results from gluino searches at the LHC Run-I have pushed the lower 
limit on its mass to the TeV scale [9-12]. While all of this indicates that SUSY lies at the 
TeV scale, such a heavy sparticle mass spectrum might spoil the naturalness of the MSSM 
by requiring excessive fine-tuning for generating the correct Higgs boson mass [13]. 

In the MSSM, the minimization of the tree-level Higgs potential leads to the following 
relation between the mass of the Z boson, m^, and the soft SUSY-breaking Higgs sector 
parameters [14]: 



-I? - . 


( 1 . 1 ) 


The last approximation in the above equation assumes tan/3 > 10, where tan/3 = 
with Vu being the vacuum expectation value (VeV) of the u-type Higgs doublet and 
that of the d-type one. ttiHu ^md are the soft SUSY-breaking masses of these two 





Higgs doublets and the parameter ji is the common mass paramater for the two Higgs 
superfields, originating in the superpotential of the MSSM. In order to avoid a large fine- 
tuning in Eq. ( 1 . 1 ), /i and ttiHu ought to lie in the ^ 100 GeV - 200 GeV range. 

The MSSM contains four neutralinos, Xi -45 which are the mass eigenstates resulting 
from the mixing of the fermion components of the Higgs superfields, known as the higgsi- 
nos with those of the gauge superfields, the gauginos The lightest of 

these neutralinos is a dark matter (DM) candidate when i?-parity is conserved. The phys¬ 
ical masses of these neutralinos are dependent on the soft SUSY-breaking gaugino mass 
parameters, Mi^2^ as well as the Higgs-higgsino mass parameter fi mentioned above. The 
unification of the soft gaugino masses at some very high scale implies that Mi and M 2 are 
of the same order as the gluino mass parameter M 3 at the SUSY-breaking scale [15-18]. 
Thus the exclusion limits on the gluino mass from the LHC, together with the requirement 
of naturalness, lead to a large splitting between the parameters Mi^2 and /i. This in turn 
implies small gaugino-higgsino mixing and, after diagonalization of the neutralino mass 
matrix, two of the physical neutralinos are gaugino-like while the other two are almost 
purely higgsinos, which are very close to each other in mass. In fact, for Mi, M 2 ^1.5 TeV 
and fi ^ 150 GeV, the mass-splitting, m^o — between the two lightest higgsino-like 
neutralinos is less than 5 GeV. At the same time, the lighter of the two charginos, is 
also a pure higgsino while the heavier, a gaugino. In such a scenario, the mass splitting 
between the lightest chargino and the lightest neutralino, m~± — is typically about 
half of m~o — m~o . 

A2 

To search for SUSY in the parameter space regions of the MSSM with nearly degenerate 
higgsinos is one of the major challenges for particle colliders. A lot of emphasis in this 
regard has been laid on the mono-jet, mono-photon or mono-Z searches at the future 
experiments [19-36]. However, owing to the very small signal rates as well as the statistical 
limitations, all these channels are expected to show only percent level yields at the 14 TeV 
LHC. Several studies [37-51] have suggested that the presence of extra leptons may help 
in improving the sensitivity for these processes. In [52, 53] such a compressed higgsino 
spectrum has been probed by tagging an ‘opposite sign - same flavor’ (OS/SF) lepton 
pair originating from the decay of a heavy neutralino. However, this method only works 
well when the higgsino mass-splitting is around 10 GeV or larger. This is for two main 
reasons. First, for smaller mass-splitting the two leptons produced are too soft to be 
tagged efficiently. Second, since these leptons are highly collinear, the signal events are 
diminished by the requirement to isolate them individually. 

Therefore, it is imperative to develop new methods for exploring regions of the MSSM 
parameter space which are consistent with the naturalness criteria [54-62] but which may 
have stayed hidden at the LHC so far. In this article, we discuss a method for probing the 
compressed higgsino spectrum in which the two highly collimated muons produced in the 
decays of X 2 identified as a single object. We explain the event selection procedure, 
specific to the kinematics of our signal process, which can be employed to reduce the 
backgrounds. Using some benchmark MSSM points consistent with such a scenario, we 
analyze the sensitivity that can be achieved at the 14 TeV LHC using our method. 

The article is organized as follows. In section 2, we briefly discuss the model parameter 
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Figure 1: Mass-splitting between the two higgsinos as a function of the gaugino mass 
parameters Mi and M 2 , ji — 150GeV and tan/3 = 30. 


configurations leading to the scenario of our interest. In section 3 we explain a tagging 
method for two soft and collimated muons. In section 4 we discuss our numerical results 
in detail. Finally, we present our conclusions in section 5. 


2 Nearly mass-degenerate higgsinos in the MSSM 


The tree-level neutralino mass matrix in the MSSM is written, in the basis 
as 


A4yo 


/Ml 0 

M 2 

—mw tan 9w cos j3 mw cos j3 
y mvi/tan 014/sin/3 —mvi/sin^S 


—mw tan 9w cos /3 mw tan 9w sin /3 \ 
mvi/ cos j3 —mw sin /3 


0 -fi 

-/i 0 y 


( 2 . 1 ) 


where mvi/ is the mass of the VF boson and 9w is the weak mixing angle. The above 
mass matrix can be diagonalized with an orthogonal real matrix as NM^oN"^ — 
diag(m^o, m^o, m^o, m^o), such that m^o < m^o < m^o < m^o. By assuming Mi, M 2 ^ 
|/i| in the above mass matrix, one obtains the approximate relation. 


Am 


m^o — m^o 

A2 ai 


rW , 6>iy 

M 2 ^ Ml 


( 2 . 2 ) 


Similarly, using also the chargino mass matrix, one gets m-± — m^o « (ignoring the 
terms proportional to l/tan/3). The neutralino mass matrix in Eq. (2.1) is subject to 
higher order corrections. The diagonalization of the mass matrix in which such corrections 
have been included (at a certain perturbative order) can be conveniently done numerically 
using publicly available SUSY mass spectrum calculators. We used the program SUSY- 
HIT vl.4 [63] to scan over Mi and M 2 , both ranging from 0.5TeV to lOTeV. For this scan 
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Am = m^o — m^o [GeV] 


Figure 2: Lifetimes of and X 2 functions of Am, which is varied by adjusting M 
(= Ml = M 2 ). The values of fi and tan/3 used are the same as in figure 1. 


we set ji = 150 GeV, tan = 30 and the input masses of the SUSY particles other than the 
electroweakinos to very high values so that they are effectively decoupled. The resultant 
values of Am are shown in figure 1. We see that for Mi^2 ^ 1.5 TeV the mass-splitting 
between the higgsino-like y? ^ is always less than 5 GeV. Also, for such large Mi^2, 
the DM direct detection facility XENONIT [64] will not be sensitive to the y? [28, 34], 
which is the lightest SUSY particle (LSP), when its mass is less than 200 GeV. 

Since X 2 almost mass-degenerate with x^, the former can live long enough to leave 
a secondary vertex in the detector. In figure 2, we show the liftime of X 2 ^ calculated with 
SUSY-HIT, as a function of Am. We see that, for Am < 1 GeV, the lifetime of X 2 
be long enough to produce displaced vertices of order 100 fim. In fact, for Am < 0.1 GeV, 
X 2 can become collider-stable, so that it leaves the detector before decaying. However, 
such a tiny higgsino mass-splitting only occurs for Mi^2 of 0(100TeV). One also sees in 
the figure that a chargino has a slightly longer lifetime than a neutralino, which is because 
m~± — m^o as noted above. A strong limit on the chargino lifetime, shown by the 

pink/shaded region in the figure, has recently been obtained by the CMS collaboration [65] 
for Am < 1 GeV. At the LHC, a mono-jet along with a displaced vertex larger than 100 jim 
might help probe the region with Am <2.5 GeV. 

A number of studies, as noted earlier, have explored the MSSM regions with Am > 
10 GeV in decays of ^ that involve two leptons in the final state. Here we will focus on 
the splitting region 2 GeV < Am < 5 GeV for the decay process 

X 2 ^ ^ Xi^^^ 5 (2-3) 

with the X 2 produced via pp X 2 X 1 + Iri our case, due to the small mass-splitting 
between the higgsinos, the two leptons are generally very soft. Therefore, we only consider 
muons in the final state on account of a much cleaner background as well as a much higher 
trigger efficiency in their case compared to those for tans or electrons. 
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Figure 3: The separation between the two muons coming from the Z*. The dotted 

vertical line shows that with the conventional cut, = 0.3, the muons in the signal 

process can not be isolated. 


In figure 3 we show the separation, Ar]^ + Ac/)^ (with r] being the pseudo¬ 

rapidity and (j) being the azimuthal angle), between the two final-state muons {pt{i^) > 
5 GeV), for Am — 3,4,5GeV. Evidently, the usual isolation criteria for a single lepton, 
Ai?max _ Q 3 ^ remove a large number of the signal events. Thus we need to use an 
unconventional reconstruction method for probing such collimated muons and establishing 
our signal over the SM backgrounds. 

3 Reconstructing collimated muons 

In order to probe the highly collinear muons produced for a very small Am, we cluster 
them together into one object, /icoi, during our simulation of the higgsino pair-production 
process. This method is similar in concept to the identification of a depton-jet’ [66-71] and 
has already been used recently in analyses of the decays of a light dark photon or of a light 
scalar or pseudoscalar (< 3 GeV) [72] into two or more soft leptons. Instead of imposing 
the conventional criterion of AR^^^ = 0.1 in order to identify the two muons coming from 
the Z* as a lepton-jet, we use a modification of the criteria described in the CMS analyses 
[72, 73] for probing collimated muons. Our method is explained below. 

• Capturing /icoP We require pT > 5 GeV for each muon in the signal, before isolation. 
In addition to this, we impose the cut m^^ < 5 GeV on the invariant mass of the 
muon pair, since we are only interested in Am < 5 GeV. 

• Isolation: To suppress the backgrounds containing muon pairs from meson decays, we 

apply an isolation criterion, /gum < 3 GeV, on pco\- The isolation parameter /gum is 
defined as the scalar sum of the transverse momenta of all additional charged tracks, 
each with pT > 0.5 GeV, within a cone centered along the momentum vector of pco\ 
and satisfying = 0.5. 
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Figure 4: Feynman graph for our signal process. The cone signifies a pair of collimated 
muons /icoi coming from the Z*. j denotes a hard QCD jet from initial state radiation. 


For the Monte Carlo simulations, we generated the parton-level signal and background 
events with MADGRAPH_aMC@NLO [74]. These events were then passed on Pythia 6 
[75] for hadronization and subsequently to the fast detector simulator Delphes 3 [76] 
interfaced with Fast Jet- v3.0.6 [77] for jet-clustering. In Delphes 3 we added a class for 
/icoi identification. The jets were clustered using the anti-fer [78] algorithm with 
set to 0.4. As a test of the implementation of our method, we first performed simulations 
for the benchmark points provided in the CMS analysis [72] and found our results to be 
within 5% of the ones presented there, in terms of signal efficiencies. 

According to [72, 73], the largest background for our signal process is the bb production, 
which has a cross section of 0(10^ pb). Although requiring 5-quarks decays, via double 
semileptonic decays, into pairs of muons (the branching ratio being around 1%) which 
are isolated, reduces this background by 0(10“^), it is still huge. We therefore apply a 
parton-level cut of pT > 200 GeV on the first leading jet in the signal and the backgrounds. 
This translates into the requirement of large missing energy, in the final state at the 
detector level, which almost entirely removes the bb background. To include QCD effects 
[79], we use the MLM-scheme to match the additional two jets [80]. We illustrate our signal 
process in figure 4. 

In figure 5 we show the combined cross section, after imposing the pT > 200 GeV cut 
on the leading jet, for our signal process, pp X 2 X 1 + two additional processes, 

pp X 2 xf^ three different values of Am. The reason for including 

the latter two processes is that the xf in the second process as well as the additional X 2 
in the third process gives very soft products which escape undetected, thus resulting in 
only $lrp in the final state and thereby mimicking the signal process. We see in the figure 
that the cross section gets considerably reduced for smaller mass-splitting. However, even 
with Am = 3 GeV for a ^ 200 GeV y?, more than 20 signal events can be obtained at the 
14TeV LHC with an integrated luminosity, £, of 3000 fb“^. 
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Figure 5: The leading order cross sections corresponding to three different values of Am, 
as a function of the x? mass. We require pT > 200 GeV for the leading jet, pT > 5 GeV for 
each muon and the detector geometry cuts to be satisfied. 


For our signal-to-background analysis, we choose three benchmark points (BPs) with 
a different value of Am each. Some specifics of these points are given in Table 1 . The 
signal cross sections given in the table have been obtained after requiring pT > 200 GeV 
for the leading jet and pT > 5 GeV for each muon. The detector geometry cuts on the 
pseudorapidities of the muons (|p| < 2.5) and the leading jet (|p| < 5.0) have also been 
imposed. 

Table 1: Some properties of the three benchmark points analyzed in this study. 



BPl 

BP2 

BP3 

m^o (GeV) 

103.0 

120.0 

150.0 

Am (GeV) 

3.0 

4.0 

5.0 

Signal cross section (fb) 

0.12 

0.165 

0.116 


3.1 Backgrounds 

After removing the bb background for mono-jet production along with soft muons and 
the main backgrounds that remain include the following. 


• V + 7 * + jets: A large results from the W ^ Iv decay or the Z ^ vv decay 
and the two collinear muons originate from the virtual photon, 7 *. To reduce this 
background, we require m^^ > 1.0 GeV and Ai?^^ >0.1. 

• rf + jets: Each tan decays into a muon and a pair of neutrinos. Due to the large boost 
in the leading jet, these muons become highly collimated. The neutrinos produced 
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are responsible for a large 

We also look into some other backgrounds, described below, which are only 0(10“^) of the 

above main backgrounds. 

• Vbb + jets and Zb + jets: This background mimics our signal when the 6 -quarks de¬ 
cay into pairs of muons via double semileptonic decays. However, it is suppressed by 
a factor ^ 10“^ after isolation. We can estimate this background from the ATLAS 
mono-jet search [81], where a set of mono-jet cuts is imposed on the backgrounds 
{pt > 280 GeV for the leading jet and > 220 GeV). The cross section for the 
Z{-^ uTy) -\- jets background is around 0.82 pb. The cross section for W -\- jets^ where 
W decays semileptonically, is around 0.6 pb. The total cross section for all these 
backgrounds thus adds up to about 1.4 pb at the 8 TeV LHC. At the 14TeV LHC, 
even if one assumes the cross section for these backgrounds to increase by a factor of 
10, it will reduce to ^ 0.3 fb after taking into account the suppression of 2 x 10“^ from 
the possibility of the 6 -jets giving collinear muons. Thus this background becomes 
much smaller than the H + 7 * + jets backgrounds, which are still around 20 fb after 
the mono-jet cuts. 

The W{^ iJiv)bb + jets background can also mimic our signal if a 6 -jet is miss-tagged 
as a muon. The 6 miss-tagging rate after passing the isolation criteria is less than 
0.005, but since this muon tends to have a large separation from a muon resulting 
from the W decay, this contribution is also small. 

• tt-\-jets: After imposing the mono-jet cut and requiring < 5 GeV, the cross 
section for this background is reduced to less than 0.1 fb. We, therefore, do not take 
it into consideration here. 

3.2 Summary of the cuts 

Below we summarize our cuts based on the discussion above. 

• Mono-jet cut: We require pt > 250 GeV for the leading jet and veto events which 
have more than three jets with pT > 30 GeV. A(/) between the leading jet and the 
second jet should be larger than 0.4. All the jets are 6 - and r-vetoed. Any events 
containing an electron with pT larger than 10 GeV are also vetoed. We additionally 
require > 250 GeV and demand exactly one pair of SF/OS muon candidates with 
each of these muons having pt > ^ GeV. 

• Basic requirements on /icoi- First we define an object pcoi as a two-muon system 

satisfying 1 GeV < < 5 GeV and 0.1 < < 0.5. Here the cut at the lower 

end of ARjjfi is to remove the main backgrounds where two muons emerge from a 
7 *. The /icoi is required to be isolated with /gum < 3 GeV. The pT of the pco\ our 
signal usually tends to be small. We therefore apply a cut pT < 20 GeV for the /icoi- 

• Cut on A(/)(/icoi, ^t)* To further remove the backgrounds, we add a cut on the Acj 
between and as 0.1 < Acj < tt/S. In figure 6 we show the A(/)(/icoi, ^r) distri¬ 
butions for the signal corresponding to our BP2 as well as the backgrounds. We note 




Unit normalized distribution 



Figure 6: The distribution of A(/)(/icoi, $t) fhe signal and the backgrounds. The signal 
corresponds to our BP2. The grey/shaded regions are the ones cut off in our kinematical 
analysis. 


that imposing the lower cut of A(/) > 0.1 removes much of the rf jets background 
and the upper cut leaves only about a third of the U + 7 * + jets background. As for 
the signal, this cut only removes less than 30% of the events. 

• Mass cut on To suit the mass of the /icoi in oui signal, we only select events 

with 1.5 GeV < < 4 GeV, cutting off also the small window, 3.0 GeV < < 

3.2 GeV, corresponding to the mass of the J/4^ resonance. 

• Kinematic cuts: Since in our signal the missing energy originates from two neutrali- 

nos, whereas the in the U + 7 * + jets background comes from a single V boson, 
we expect different transverse mass, Mt, distributions of the signal and the back¬ 
grounds [82]. In addition, we also impose a cut on ^j^/pt(ju^coi) because the cut on 
the pt of the leading jet also results in a boosted pcoi- We can use fhis cut 

instead of the leading jet pT due to the comparatively much smaller Pt(/^coi)- In 
figure 7 we show the distributions of these two variables for the BP2 signal (left) and 
the U + 7 * jets background (right). We find that an upper cut of Mt < 50 GeV 
and a lower cut of ^x/PT{fJ^co\) > 20 suppresses the background while allowing most 
of our signal events. 


4 Results of the signal-to-background analysis 

In Table 2 we show the cut-flow for our numerical simulations of the three BPs. We note 
here that while generating our signal process for each BP we also required Pt{p) > 4 GeV 
and 1 GeV < < 5 GeV at the parton level. The table shows that the total background 

cross section after the mono-jet cut is around 23 fb, which is three orders of magnitude 
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Figure 7: Mt vs ^7-/]9T(/icoi) foi* the signal corresponding to our BP2 (left) and for the 
Z 7 * background (right), after applying the mono-jet cut as well as the basic cuts on /icoi- 
We mark the cuts on Mt and ^T/Pr(Mcoi) with solid magenta lines. Thus the events 
allowed after the cuts are the ones located in the boxes traced by the lines in the upper-left 
corners of the figures. The heat map corresponds to the number of events. 


larger than our signal. We observe that the F + 7 * + jets background is much larger in 
our case than in [53]. The reason for the strong suppression of this background in [53] 
is that the two muons coming from the 7 * are individually isolated there. However, in 
our case such an isolation condition will conversely suppress the signal process due to a 
comparatively smaller mass-splitting between the two higgsinos. 

According to Table 2 , after applying all the cuts, the H + 7 * + jets background is 
clearly the largest one. The jrf background only contributes about 10% to the total since, 
as noted earlier, the lower cut on A(/)(/icoi, $t) reduces this background by nearly an order 
of magnitude. 

Table 2: Cut-flow for our signal and background processes. The cross sections are in fb 
and S/B and S/\fB given in the last two rows correspond to C = 3000 fb“^ at the 14TeV 
LHC. 


Cuts 

WTj 

ZTj 

jrr 

Total BKG 

BPl 

BP2 

BP3 

Mono-jet 

8.057 

8.82 

6.674 

23.0 

0.052 

0.072 

0.056 

Basic /icoi 

0.753 

1.05 

0.314 

2.1 

0.041 

0.042 

0.028 

^0(/^COl5 $t) 

0.288 

0.324 

0.035 

0.65 

0.028 

0.030 

0.020 


0.106 

0.118 

0.024 

0.248 

0.017 

0.023 

0.015 


0.037 

0.044 

0.011 

0.092 

0.013 

0.016 

0.010 

S/B 

0.14 

0.17 

0.11 

S/^/B (a) 

2.4 

2.9 

1.85 


We also note in the table that the highest significance we obtained for the 14 TeV LHC 
with C — 3000 fb“^ is ^ 3 cr and corresponds to our BP2, while S/B for this point is 17%. 
We point out here that although both BPl and BP3 give signal cross sections similar to 
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the one obtained for BP2, the obtained significance is smallest for BPS. The reason is 
that for BPS Am = 5GeV, so that i?t(Mcoi) tends to be a little larger than our chosen 
strong upper cut on pr(Mcoi)* Another factor is the < 4GeV cut, which also removes 

some signal events in the case of BPS, but not in the case of BP2. Although we can relax 
the upper cut on Pt(Mcoi) fo around 25GeV, it will also result in larger backgrounds and 
thus the statistical significance for BPS will not improve. We should also point out here 
that for BPl, where Am = SGeV, since m^^^^ should be less than S GeV, changing the 
upper cut on m^^^^ from 4 GeV to S GeV would enhance our signal significance. However, 
in the experimental searches the true value of the higgsino mass-splitting is unknown. We 
therefore retain the upper cut of 4 GeV, which suits most of the Am range that we are 
interested in. 

Finally, using our method the minimum Am we have managed to explore is ^ SGeV. 
This is because of the requirement of pr > 5 GeV for the muons. By imposing a lower cut, 
Pt > 4 GeV [8S], for the muons, a good sensitivity to even smaller values of Am can be 
achieved. 

5 Conclusions 

In this article, we have analyzed the possibility of probing natural SUSY scenarios with 
a highly compressed higgsino mass spectrum at the 14TeV LHC, using the collinearity 
between the two muons produced in such scenarios. We have found that a statistical 
significance of up to 3 a as well as S/B up to 17% can be obtained for Am = m^o — m^o = 
4 GeV and a ^ 120 GeV Xi with an integrated luminosity of 3000 fb“^ at the LHC. In fact, 
by using our analysis method but further lowering the cut we imposed on the pT of muons, 
MSSM parameter space regions with Am even lower than 3 GeV can be explored. 
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